Scalar Casimir effect for D-dimensional spherically symmetric 

Robin boundaries 



o 
o 

(N 



(N 
> 

in 

oo 



X 



Aram A. Saharian* 
February 1, 2008 



Department of Physics, Yerevan State University, 
1 Alex Manoogian Street, 375049 Yerevan, Armenia 



The vacuum expectation values for the energy-momentum tensor of a massive scalar field 
with general curvature coupling and obeying the Robin boundary condition on spherically sym- 
metric boundaries in D-dimensional space are investigated. The expressions are derived for the 
regularized vacuum energy density and radial and azimuthal stress components (i) inside and 
outside a single spherical surface and (ii) in the intermediate region between two concentric 
spheres. A regularization procedure is carried out by making use of the generalized Abel-Plana 
formula for the series over zeros of cylinder functions. The asymptotic behavior of the vacuum 
densities near the sphere and at large distances is investigated. A decomposition of the Casimir 
energy into volumic and surface parts is provided for both cases (i) and (ii). We show that 
the mode sum energy, evaluated as a sum of the zero-point energy for each normal mode of 
Q_i| frequency, and the volume integral of the energy density in general are different, and argue that 

this difference is due to the existence of an additional surface energy contribution. 
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1 Introduction 

Historically the investigation of the Casimir effect (for a general introduction, see [l], ^, |3|) 
for a perfectly conducting spherical shell was motivated by the Casimir semiclassical model 
of an electron. In this model Casimir suggested that Poincare stress, to stabilize the charged 
particle, could arise from vacuum quantum fluctuations and the fine structure constant can 
be determined by a balance between the Casimir force (assumed attractive) and the Coulomb 
repulsion. However, as has been shown by Boyer the Casimir energy for the sphere is 
positive, implying a repulsive force. This result was later reconsidered by a number of authors 
[Hi Hi @]- More recently new methods have been developed for this problem including direct 
mode summation techniques based on the zeta function regularization scheme ]|]-f2(| (for similar 
considerations in the case of a dielectric ball see, for instance, the references given in pi]]). 

Investigation of the dimensional dependence of physical quantities in the Casimir effect is 
of considerable interest. In particular, for periodic boundary conditions this is motivated by 
the idea of using the Casimir effect as a source for dynamical compactification of the extra 



dimensions in Kaluza-Klein models. In [22] the Casimir energy is derived in a general hyper- 



cuboidal region and for various types of boundary condition and field. For a spherical shell the 
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Casimir effect in an arbitrary number of dimensions is analyzed in || [T(| for a massless scalar 
field satisfying Dirichlet and a special type of Robin (corresponding to the electromagnetic TM 



modes) boundary conditions using the Green's function method (see also p3[ ) and in [19] for 
the electromagnetic field and massless scalar and spinor fields (for the 4D fermionic Casimir 
effect see |2j, |25|, p6| ) with various boundary conditions on the basis of the zeta regularization 
technique. 

However, the most of the previous studies on spherical geometry were focused on global 
quantities such as the total energy and the force acting on a shell. Investigation of the energy 
distribution inside a perfectly reflecting spherical shell was made in [27] in the case of QED and 
in [28] for QCD. The distribution of the other components for the energy-momentum tensor of 
the electromagnetic field inside as well as outside the shell can be obtained from the results of 
[ p^ , 30]. In these papers the consideration was carried out in terms of the Green's function (for 
this formalism see f7|). Another efficient method for extracting finite parts from the vacuum 
expectation values of the local physical observables in the case of plane boundaries is based 
on the Abel-Plana summation formula (see, e.g., |IJ). We have generalized this formula to 



include curved boundaries [31], g^] . In |33|, £4| (see also (3^] ) the calculations of the regularized 
vacuum expectation values for the electromagnetic energy-momentum tensor inside and outside 
a perfectly conducting spherical shell and in the region between two concentric spheres are based 
on the generalized Abel-Plana formula. 

In this paper the vacuum expectation values of the energy-momentum tensor are investigated 
for a massive scalar field with general curvature coupling parameter £, satisfying the Robin 
boundary condition on spherically symmetric boundaries in D-dimensional space. As special 
cases they include the results for the Diriclet, Neumann, TM, and conformally invariant Hawking 
boundary conditions. Robin type conditions also appear in considerations of the vacuum effects 
for a confined charged scalar field in external fields (see, for instance, |35|]) and in quantum gravity 
[|3q] . In addition to describing the physical structure of the quantum field at a given point, the 
energy- momentum tensor acts as the source of gravity in the Einstein equations. It therefore 
plays an important role in modeling a self-consistent dynamics involving the gravitational field 
[j37fl. To the author's knowledge all previous investigations on the spherical Casimir effect for 
scalar fields were concerned mainly with global properties (interior and exterior energies, vacuum 
stress on a sphere). The massless scalar field with Dirichlet boundary condition is considered 
in |}8], 11, [l2|, H|, [l6|, [l9| and for Neumann and Robin boundary conditions in 10, 16. |X9fl , 



The case of the massive scalar field is investigated in [Q5| (for the heat-kernel coefficients and 
determinants, see 14] and references therein). In |39|, p0| asymptotic expansions for the 
renormalized energy-momentum tensor are developed near an arbitrary smooth boundary in the 
case of conformally and minimally coupled 4D massless scalars with Dirichlet, Neumann, and 
Robin boundary conditions. Here we consider the scalar vacuum inside and outside a single 
spherical shell, and in the intermediate region between two concentric spherical shells. For 
the latter geometry the general case is investigated when the constants in the Robin boundary 
condition are different for inner and outer spheres. To evaluate the corresponding field products 
we use the mode sum method in combination with the summation formulas from Ref. |3l|] 
(see also |32|). For scalars with general curvature coupling the essential point is the relation 
between the mode sum energy, evaluated as a renormalized sum of the zero-point energies for 
each normal mode of frequency, and the volume integral of the renormalised energy density. For 
flat spacetime backgrounds the first quantity does not depend on £, since the normal modes 
are the same for fields with different values of this parameter (for instance, for minimal and 
conformal couplings). Nevertheless, the corresponding energy-momentum tensor depends on £ 
and, in general, this is the case for the vacuum energy distribution and hence for the integrated 
vacuum energy. As a result, as was mentioned in [39], the mode sum energy and integrated 
energy in general are different. (Note that in most of the papers referred to above the first 
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quantity is considered.) Below for the geometries under consideration we calculate both these 
quantities and argue that this difference is due to the existence of an additional surface energy 
contribution to the total vacuum energy, and the Casimir energy decomposition into volume 
and surface parts is provided (for a similar consideration in the case of parallel plate geometry 
see ||). 

We have organized the paper as follows. In the next section we consider the vacuum inside 
a sphere and derive formulas for expectation values of the energy density and stresses. Section 
H is devoted to the corresponding global quantities, such as the total interior Casimir energy 
and vacuum force acting on a sphere. We show that the interior vacuum energy contains two 
parts: volume and surface ones. The vacuum expectation values of the energy- momentum tensor 
(EMT) for the region outside a sphere are considered in section f|. The expressions for the total 
energy and force acting on a sphere from outside are derived. The total outside Casimir energy 
is decomposed into surface and volume parts. Further in this section a spherical shell with zero 
thickness is considered and the corresponding global quantities including interior and exterior 
parts are evaluated. In section [| we consider the vacuum energy-momentum tensor between 
two concentric spheres and section ||] is devoted to the global quantities in this region. Section 
[j] concludes the main results of the present paper. In the appendix we derive the summation 
formula over zeros of the Bessel functions combination on the basis of the generalized Abel-Plana 
formula. The vacuum expectation values in the region between two spheres contain this type of 
series. 



2 Vacuum EMT inside a spherical shell 

Consider a real scalar field ip with curvature coupling parameter £ in D-dimensional space 
satisfying the Robin boundary condition 



{Ax + Bm*Vi) tp(x) = 



(2.1) 



on the sphere, where A\ and B\ are constants, n % is the unit inward normal to the sphere, 
and Vj is the covariant derivative operator. Of course, all results in the following will depend 
on the ratio A\/B\ only. However, to keep the transition to the Dirichlet and Neumann cases 
transparent we will use the form (2.1). The corresponding field equation has the form 

1 



(ViV 1 + m 2 + £R) ip = 0, ViV J 



-g 



-A 



99 ik d k ), 



(2.2) 



where R is the scalar curvature for the background spacetime. The values £ = 0, and £ = £ c 
with £ c = (D — 1)/4D correspond to the minimal and conformal couplings, respectively. Here 
we will consider the case of flat spacetime. The corresponding metric energy-momentum tensor 
is defined as (see, e.g., f37j) 



T ik = (1 - 20di<pd k ip + (2£ - l/2)g lk d l ^dw ~ 2£^V 4 V^ + (1/2 - 2£)m 2 <^ 2 . (2.3) 
It can be seen that by using the field equation this expression can also be presented in the form 



T ik = di(pd k <p + 



and the corresponding trace is equal to 



(2.4) 



(2.5) 
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By virtue of Eq.(2.4) for the vacuum expectation values (VEV's) of the EMT we have 



{0\T lk (x)\0) = lim ditf k (0\<p(x)ip(x')\0) + 



(0|v9 2 (x)|0), (2.6) 



where |0) is the amplitude for the vacuum state. Note that the VEV (0\<p(x)ip(x')\0) = G + (x, x') 
is known as a positive frequency Wightman function. In Eq. (|2.6|) instead of this function one 
can choose any other bilinear function of fields such as the Hadamard function, Feynman's 
Green function, etc. The regularized vacuum EMT does not depend on the specific choice. The 
expectation values (|2.6|) are divergent. They are divergent in unbounded Minkowski spacetime 



as well. In a flat spacetime the regularization is performed by subtracting from Eq.(2.6) the 
corresponding Minkowskian part: 



{Tik(x)) reg = (0|T ifc (x)|0) - {0 M \Tik(x)\0 M ) = 9 ik {<p(x)<p(x')) 



reg - 



(2.7) 



where |0m) denotes the amplitude for the Minkowski vacuum and the form of the operator 9i k 
directly follows from Eq,([2.6j). Therefore the finite difference between two divergent terms in 
Eq.(2.7) can be obtained from the corresponding difference between the Wightman functions, 
by applying a certain second-order differential operator and taking the coincidence limit. To 
derive the expression for the regularized VEV of the field bilinear product we will use the 
mode summation method. By expanding the field operator over eigenfunctions and using the 
commutation rules one can see that 

(0\<p(xMx')\Q) = £ <Pa(z)<P*atf), (2-8) 

a 

where {cp a (x), (p* a (x')} is a complete set of positive and negative frequency solutions to the field 



equation ( |2.2j ), satisfying boundary condition (2.1). 

In accordance with the symmetry of the problem under consideration we shall use hyper- 
spherical polar coordinates (r, 1?, (f>) = (r, 0%, 62, ■ ■ ■ G n , (ft), n = D — 2, related to the rectangular 
ones (xi, X2, xd) by (see, for instance, Section 11.1) 



x\ 

XD-l 



r cos 0i, , X2 = r sin 0i cos #2, ••• 
r sin 0i sin 02-" sin 9 n cos <f>, xjj 



x n = r sin 9\ sin 02... sin n _i cos 9 n (2.9) 
r sin 0i sin 02 • • • sin 9 n sin <p, 

where < 9 k < tt, k = 1, n and < 4> < 2n. In the hyper spherical coordinates for the region 
inside the sphere the complete set of solutions to Eq^^), regular at the origin, has the form 

-iu)t 



ip a {x) = P a r- n/2 J u (r^uj 2 - m 2 )Y(m k ;tf, 0)e~ tuJt , v = l + n/2, I = 0, 1, 2, 
where m/% = (mo = I, mi, ■ ■ ■ , m n ), and mi,m^, ■ ■ ■ , m n are integers such that 
< m n _i < m n _2 < • • • < mi < I, -m n -i < m n < m n -\, 



(2.10) 



(2.11) 



J u {z) is the Bessel function, and Y{m^ \ ??, 4>) is the surface harmonic of degree I (see |42j, Section 
11.2). This last can be expressed through the Gegenbauer or ultraspherical polynomial Cp(x) 
of degree p and order q as 



k=l 



The corresponding normalization integral is in the form 



\Y(m k ;$,<P)\ 2 dn = N(m k ). 



(2.12) 



(2.13) 
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The explicit form for N(rrik) is given in |^] , Section 11.3, and will not be necessary for the 
following considerations in this paper. From the addition theorem [42], Section 11.4, one has 

1 



^Y {mM) Y*^>A>) = ^^ 2 (cos ( 



(2.14) 



where S D = 2vr D/2 /Y{D/2) is the total area of the surface of the unit sphere in D-dimensional 
space, 9 is the angle between directions ($,</>) and and sum is taken over the integer 

values m k , k = 1, 2, . . . , n in accordance with Eq.( |2.11| ). 

The coefficients (3 a in Eq.( |2.10| ) can be found from the normalization condition 



\(f a (x)\ dV 



1 

2u' 



(2.15) 



where the integration goes over the region inside the sphere. Substituting eigenfunctions (2.10), 
and using the relation ( |2.13| ) for the spherical harmonics and the value for the standard integral 
involving the square of the Bessel function, one finds 



0l 



A 



N(irik)um 



T u (Xa) 



with the notations 



uJ- 



Tu(z) 



(z 2 - ^)J 2 (z) + z 2 J' 2 (z)' 



(2.16) 



(2.17) 



From boundary condition (2^1) on the sphere surface for eigenfunctions (|2.10| ) one sees that 
the possible values for the frequency have to be solutions to the following equation 



AJ v (z)+BzJ' v (z) =0, z = Xa, A = A 1 + B 1 n/2a, B = -B 1 /a. 



(2.18) 



It is well known (see, e.g., [f43| , |44|| ) that for real A, B and v > — 1 all roots of this equation 
are simple and real, except the case A/B < —v when there are two purely imaginary zeros. 
Let us denote by X U)k , k = 1,2, ... , the zeros of the function AJ v {z) + BzJ' v {z) in the right 
half plane, arranged in ascending order of the real part. Note that for the Neumann boundary 
condition A/B = 1 — D/2, for the TE and TM electromagnetic modes B = (Dirichlet) and 
A/B = D/2 — 1, respectively, and for the conformally invariant Hawking boundary condition 
@ one has A/B = 1/D- D/2. 

Substituting Eq.( |2.10| ) into Eq. Q2.8j ) and using addition formula ( 2.14| ) for the spherical har- 
monics, one obtains 



{0\<p(x)<p&)\0) 



( rr l\-n/2 00 



naSo 



(2.19) 



1=0 



; '= ' = J u (\ Ujk r/a)J u (X Utk r /a)e V 



k=l 



To sum over k we will use the generalized Abel-Plana summation formula [31 

Yu(z) 



2j2 T A^,k)f{K,k) = / f(x)dx+ -Res, =0 /(z 
k=i ^° 



Ju(z) 



Ux) 



/Ore™ /2 ) + e w ' 



f(xe 



(2.20) 
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where, following pl| , for a given function F(z) we use the notation 

F(z) = AF(z) +BzF'(z). 



(2.21) 



Formula ( 2.20| ) is valid for functions f(z) analytic in the right half plane and satisfying the 
conditions 



2uiri 



f(z)+o(z 2 ^ 1 ), 



\f{z)\<e{x)e c \y\ c<2, s(x) - 0, 



z - 

X 



0, 



oo. 



(2.22) 
(2.23) 



This formula can be generalized in the case of the existence of purely imaginary zeros for the 
function J v {z) by adding the corresponding residue term and taking the principal value of the 
integral on the right (see [j32|). However, in the following we will assume values of A/B for 
which all these zeros are real. 

As the function f(z) in Eq.( 2.20 ) let us choose 



/(*) 



\/ z 2 + m 2 a 2 



J u (zr/a)J u (zr' '/ a ) e V z2 /« 2 +m 2 (t'-t). 



(2.24) 



This function has branch points on the imaginary axis. As has been shown in (3l], |32| the formula 
fl2.20p can be used for functions having this type of branch point as well. As f(z) ~ z 2u+1 , z — > 
condition ( 2.22 ) is satisfied. Using the asymptotic formulae for the Bessel function (see, e.g., 
p4[| ) it is easy to see that condition (2.23) is satisfied if r + r' + \t — t'\ < 2a. (Note that this 
condition is satisfied in the coincidence limit r = r', t = t' for interior points, r < a.) Assuming 
that this is the case and applying to the sum over k in Eq.(HJ|) formula ( gjg ) with f(z) from 
Eq.( |2.24] ) one obtains 



(0| v (xMx')|0) 



E 



21 + n 



C, n/2 (cos0) 



2naS D ^ {rr'yl 2 1 



f(z)dz 



(2.25) 



o 



00 ^ K v (z)I u (zr/a)I u (zr'/a) , 
az z = — cosh 

ma Iu(Z) Vz 2 



m 2 a 2 



^z 2 /a 2 -m 2 (t' - 1) 



where we have used the result that the difference of the radicals is nonzero above the branch 
point only and introduced the modified Bessel functions. The contribution of the term in the 
first integral to the VEV does not depend on a, whereas the contribution of the second one tends 
to zero as a — > oo. It follows from here that the first term is the corresponding function for the 
unbounded Minkowski space: 



(0 m \(P(x)<p(x')\0m) 



1 



E 



21 + n 



Cr /2 (cos#) 



2nS D (rr'W 2 1 



ze iVz 2 +m' 2 (t'-t) 

dz , — J y {zr)J v {zr'). 



V z 2 + m 2 



(2.26) 

This can be seen also by direct evaluation. Indeed, the sum over I can be summarized using the 
Gegenbauer addition theorem for the Bessel function J23] . The value of the remaining integral 
involving the Bessel function can be found, e.g., in |45| ]. As a result we obtain the standard 
expression for the D-dimensional Minkowskian Wightman function: 



(0 m \v(x)(p(x')\0m) 



in 



(D-l)/2 



K{D-i)/2 (W(x-x') 2 -(£-*') 2 



l£ 



( 2vr )(D+l)/2 [(x _ x/)2 _ (t _ t ,y + - £ ] 



(2.27) 



where x = (t, x), e > for t > t', and e < for t < t'. Using the Feynman Green function in 
Eq.(|2.6[) we would obtain the same expression (2.27) with e > 0. 
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Hence we see that the application of the generalized Abel-Plana formula allows us to extract 
from the bilinear field product the contribution due to the unbounded Minkowski spacetime. 
Combining Eqs. fl2.19 ) and (2.26) for the regularized Wightmann function one obtains 



((p(x)ip(x')) 



V r,y 



(0\(p(x)<p(a/)\0) - (0m\<p(xMx')\0 m ) 



1 ^ 21 + n n/2 . 



00 K v (za) I u {zr)I v (zr>) 
az z^et-, — ; - cosh 



I v {za) yfz 2 



mrS D (rr') n / 2 1 
Vz 2 - m 2 (t' - t)] . 



(2.28) 



by this function 46]. Performing the limit x' — > x and using the value 



Let us recall that the response of the particle detector at a given state of motion is determined 

(2.29) 



a 



n/2. 



1) 



T(l + n) 



T{n)l\ 



(see 1 42], Section 11.1) for the regularized field square we get 

K v {z) zl 2 (zr/a) 

u x I az 

/-o 

where 



<^ 2 (r)> 



rcg 



nar n S D f-' J ma I v 



(z) \J z 2 — urn 



2 n 2 



Di = (21 + D - 2) 



r(z + D - 2) 



(2.30) 



(2.31) 



Y(D - 

is the degeneracy of each angular mode with given I. 

The VEV for the EMT can be evaluated by substituting Eqs.(^28|) and ( ^30| ) into Eq.(^7|), 
where the operator 9^ is defined in accordance with Eq.fl2.6l). From the symmetry of the 
problem under consideration it follows that (Tjfc) reg is a combination of the second rank tensors 
constructed from the metric g^, the unit vector t l in the time direction, and the unit vector n % 
in the radial direction. We will present this combination in the form 



(T, 



i )reg 



eiii k +priin k + p± (iii k - rnn k 



diag(e, -p, -p±, . . . ,-p±) , 



(2.32) 



where the vacuum energy density e and the effective pressures in radial, p, and azimuthal, 
p±, directions are functions of the radial coordinate only. As a consequence of the continuity 
0, these functions are related by the equation 



equation Vfc(Zf) reg 



r£ + (D-l){p-p ± ) 



0. 



(2.33) 



From Eqs.fl2~7l), fl2.28|) , and (I2.30D for the EMT components one obtains 

z K v {z)F^[I u (zr/a)] 



2vra 3 



dz z 6 



I u {z) \/z 2 - m 2 a 2 



e, p, p ± , r < a, (2.34) 



where for a given function f(y) we have introduced the notations 



F^ £) [f(y)] 



(i-46 



f 2 {y)--f(y)f\y) + 



y 



v 1 1 + 4£ - 2(mr/yY 



1-4^ 



f(y) 



(2.35) 



4 p Hf(y)} = f' 2 (y) + -f(y)f'(y) 

y 

4 p±) lf(y)] = (^-i)f' 2 (y)-jf(y)fUy) + 



i + 



j/2 i t n' 

- ] f(y), £i=4(n + l)£-n (2.36) 



4£-l + 



^ 2 (l+gi) + gi| 
(re + l)y 2 



f(y)- (2-37) 
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It can easily be seen that components ( |2.34|) satisfy Eq.(2.33) and are finite for r < a. The 
formulas ( p.34 ) may be derived in another, equivalent way, introducing into the divergent mode 



sum 



(0\T ik {<p(x),(p(x)}\0) = ^2T ik {(p a (x),ip* a (x)} 



(2.38) 



a cutoff function and applying the summation formula ( 2.20D . The latter allows one to extract 
the corresponding Minkowskian part by a manifestly cutoff independent method. 



At the sphere center the nonzero contribution to VEV (2.34) comes from the summands with 
I = and / = 1 and one has 



e(a,0) 



1 



2^/2+ir(ZV2) Jm v 7 ^ 



00 z D+1 dz 



V z J Id /2-i(az) loniaz) 



(2.39) 



p(a,0) = p±{a,0) 



2 D ir D / 2 + 1 DT(D/2) 



00 z D+l dz 



tt ^ K D/2~i(az) K D/2 (az) 
l€l - 2 ) — 7~T + ?1 



Id/2-i{o>z) 

Note that for the conformally coupled massless scalar 

e(o,0) = Dp(a,0). 



lD/2(az) 



(2.40) 



(2.41) 



This can also be obtained directly from the zero trace condition. The results of the corresponding 
numerical evaluation for the Dirichlet and Neumann minimally and conformally coupled scalars 
in D = 3 are presented in Fig. In the limit of large mass, ma S> 1, using the asymptotic 










/ c \ 






b 




-0 






-0 






-0 












-0 


d / 










-0 




ma 





0.01 



-0.01 



-0.02 



-0 . 03 



-0.05 




1 2 3 4 5 



Figure 1: The Casimir energy density, a D+1 e, (a, c) and vacuum pressure, a D+1 p, (b, d) at 
sphere center for minimally (left) and conformally (right) coupled Dirichlet (a, b) and Neumann 
(c, d) scalars in D = 3 versus ma. 



8 



formulas for the modified Bessel functions for large arguments and the value of the integral 



s/z 2 - n 2 



dz 



.^m^m+i^i^) 



0(2//)' 



(2.42) 



from Eqs. (|3D and ( fOol) one obtains 
e(a,0) 



L> , , (4£ - l)(am) D+1 / 2 e- 2am , r 



D — 1 



2 D 7r( D - 1 )/ 2 T(D/2)a D + 1 



(2.43) 



The expectation values (2.34) diverge at the sphere surface, r — > a (note that for a given Z 
the integrals over z diverge as (a — r) -2 ). The corresponding asymptotic behavior can be found 
using the uniform asymptotic expansions for the modified Bessel functions, and the leading 
terms have the form 



V 



e ~ 



(D-i)r((g + i)/2)(e-e c 

2 D 7r (D+l)/2 a ( a _ r )D 

z>r((D + i)/2)(£-e 



1 



2 ^p+i)/2( a _ r ) 



1 



2feo) , 

25bo) • 



(2.44) 
(2.45) 



These terms do not depend on mass or Robin coefficient Bi, and have opposite signs for Dirichlet 
and Neumann boundary conditions. Surface divergences in renormalized expectation values for 
EMT are well known in quantum field theory with boundaries || 39, S(J. In the case of D = 3 
massless fields the corresponding asymptotic series near an arbitrary smooth boundary are 
presented in [40[ |. In particular, for D = 3, £ = from Eqs. fl2.44 ) and ( |2.45 ) we obtain the 
leading terms given in f39fl for the minimally coupled Dirichlet and Neumann scalar fields. Taking 
the limit a, r — * oo, a — r = const, one obtains the leading terms for the asymptotic behavior 
near the single plate (see, for instance, |j41| ). For a conformally coupled scalar the coefficients for 
the leading terms are zero and e,p± ~ (a — r)~ D , p ~ (a — r) 1_D . In general asymptotic series 
can be developed in powers of the distance from the boundary. The corresponding subleading 
coefficients will depend on the mass, Robin coefficient, and sphere radius. 

In the case D = 1 we obtain VEV's for the one - dimensional segment —a < x < a. Due 
to the gamma function in the denominator of the expression for D t in Eq.(g3g), now the only 
nonzero coefficients are Dq = D\ = 1. Using the standard expressions for the functions ^ = bi/2(- 2; ) 

and K± 1 / 2 (z) we can easily see that F±y 2 [l±i/ 2 (z)] = ±2/{irz), and hence the vacuum stresses 
are uniform: 



P 



7T 



z dz 



m z e 4az ( Bxz-Ax 



(2.46) 



1 



For B\ = and A\ = we obtain from here the standard results for Dirichlet and Neumann 
scalars. Unlike the stress distribution the energy distribution is nonuniform. The corresponding 
expression directly follows from Eq. (|2.34 ): 



dz- 



4az ( Bxz-Aj 



B 1 z+A 1 



1+1 + 



Biz 



4^-1 
1 -m 2 /z 2 J B 1 z + A 1 



-e 2az cosh2zr 



(2.47) 



Note that for the case of massless field the expressions ( p. 46 ) and ( |2.47| ) coincide with the general 
formulas derived in |41] for two - plate geometry if one takes there (3\ = = B\jA\ and D = 1. 
In D = 1 our assumption that all zeros for J u (z) are real corresponds to B\jA\ < 0. 
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3 Total Casimir energy and vacuum forces inside a sphere 

In the previous section we considered the scalar vacuum densities and stresses inside a sphere. 
Here we will concentrate on the corresponding global quantities. First of all note that by 
using formulae and (U) the unregularized 00-component for the vacuum EMT may be 

presented in the form 



oo oo 



/>■ — j \iKk + m2a2 



f!) 0) [Ju{K, k r/a)} , 



(3.1) 



where we have introduced the notation 



f? } [f(y)\ = (i - 40 



f' 2 (y)--f(y)f'(y) + [^ 



n 



y 



v 2 4£ + 1 



y 2 4^-1 



f(y) 



mr 



+ -MyJ f(y). (3-2) 



By using the standard integrals involving the Bessel functions for the total volume energy inside 
a spherical shell one finds 



2a 



_. oo oo 



m 2 a 2 



1=0 k=l 



1 + 



(4£ - l)^iBA 2 



v.k 



(X 2 +m 2 a 2 )[A 2 + B 2 (X 2 u , - u 2 )} 



(3.3) 



As we see in the general Robin case this energy differs from the total vacuum energy inside a 
sphere: 

^ oo oo 

Eia = 2 E A 5^w v>fcl u v>k = ^l k /a 2 + m 2 . (3.4) 

Z=0 k=l 

We will interpret this difference as a result of the presence of an additional surface energy 
contribution: 



E: 



(surf) 



(4£ - \)A X B 



2a 



oo oo 



A 



such that 



z=o 



^in - &- m + • 



(3.5) 



(3.6) 



The surface energy (|3.5| ) can be obtained independently by integrating the corresponding surface 
energy density. To see this let us note that there is a surface energy density contribution to the 
total energy density in the form |4C] 



r (surf) 

[ oo 



-(2£ - l/2)5(x,dM)ipn i d i ip, 



(3.7) 



where 5(x,dM) - is a "one sided" ^-distribution. By using Eq.( |2.19| ) the corresponding VEV 
takes the form 



(0|T ( r f) (x)|0) = (2t-l/2)5(x,dM)(d r (0\v(rMr')\0)) r ,__ 



S(r-a + 0) 
2aSijr 



D-l 



oo oo 



K, t TA\t) Mz) \n Jij(z) _ zJ , (z) j ^ 2 = A r/o(3 . 8) 



A 2 , i. + m 2 a 2 



;=o fe=i y A„ )fc 

Integrating this density over the region inside a sphere, we obtain the corresponding surface 
energy: 

Et f) = fd D x(0\T^\x)\0). (3.9) 
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It can easily be seen that this expression with Eq,(|3.8;) coincides with Eq,( |3.5| ), 

The subtracted surface energy density can be obtained by applying to the sum over k in 



Eq.( |3.8|) the summation formula (2.20) and omitting the term coming from the first integral 
on the right of this formula. An alternative way is to use in Eq.(|3.8[) the subtracted product 



{ty{r)(p{r')) SUB from Eq. fl2.28Q instead of (Q\<p(r)<p(r')\0) . As a result one obtains 



T, 



(surf) 



00 



(X) 



SUB 



5(r - a + 0) 



l-4£ 



oo „ c 
1=0 Jm 



zdz 



K v {z) 



->2 n 2 



Iu{z) 



F s [I v (zr/a)}, (3.10) 



where for a given function f(y) we have introduced thenotation 



F s [f(y)]=f(y) yf'(y)--f(y) 



Integration of this energy density gives the total subtracted surface energy: 



E, 



(surf) _ 1 - 4g 

2na 



1=0 



zdz 



ma V ' Z 2 



i 2 ri 2 



K " {Z) FMZ)\. 



(3.11) 



(3.12) 



Integrating the energy density q = e in Eq. fl2.34j ) over the volume inside a sphere we obtain 
the corresponding subtracted volume energy. Using the result that for any modified cylinder 
function Z u (y) = c\I v {y) + c<iK v {y) we have the formula 



J drrF^[Z v {zr)\ = -^F v [Z u (zr)} 



with 



F v [f(y)] = (i-^)F s [f(y)] + (y 2 



one obtains 



E; 



(vol) 



1 °° 



27TO 



1=0 



m a 



KJz] 



zdz 



Iv(z) y/z 2 - m 2 d- 



(3.13) 
(3.14) 
(3.15) 



Now, using relations ( |3.12 ) and ( |3.15| ) for the total Casimir energy inside a sphere, we get 

4 = ~EA r-Vz 2 -m 2 a 2 ^44F t [I u (z)], (3.10! 

27711 T~?> -In;, '- 



1=0 



h{z) 



where for future convenience we have defined 

Ft[f(y)] = y 2 f' 2 (y)-{y 2 + v 2 )fl(y). 



(3.17) 



As we see the dependences on the coupling parameter £ in volume and surface energies canceled 
each other and the total energy does not depend on this parameter. We might expect this result 
as the eigenmodes do not depend on £. With the help of the recurrence relations for l v {z) it 
can be seen that for the case D = 3 the Dirichlet massless scalar field interior energy (3.16) 
coincides with the result given in [38]. The force acting per unit area of the sphere from inside 
is determined as 



2ira 



D 



1 00 r 
i=o Jn 



z 3 dz 



KJz) 



1 1 1 a V ' Z 2 



m 2 a 2 I v {z) 



with the notation ( |2.36| ). Note that one has the relation 

z 2 F^=F t [f{z)]+i l F s [f{z)\. 



(3.18) 



(3.19) 
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For the massless scalar this leads to the formula 



a D S D F in = + - £ c )4T f) . (3.20) 



This relation can also be derived directly from continuity equation ( 2.33| ), taking into account 
the formula 

(D-l)p ± = £ -p -#(£-&) V;VV)re g , (3.21) 

where we have used expression (|2.5| ) for the EMT trace in the case of a massless field. Sub- 
stituting Eq,( |3.2l|) into Eq.( 2.33| ) and integrating over the interior region we obtain Eq.(3.20). 



Note that for the minimally coupled scalar (£ = 0) this formula may be presented in the form 
of the standard thermodynamic relation 

d ^(voi) = _ F . ndvD + a . mdsD} (3.22) 

where vr> = Sdo. d /D and sd = Sua 0-1 are the volume and surface area for a L>-dimensional 
ball with radius a, and a- m = E^*^ /sp is the inner surface tension. 



The expressions (|3~T2|) , (PB , (gig) and (gig) are divergent in the given forms. Moreover, 



in the sums over I the individual integrals diverge. To extract finite results we can apply 
the procedure already used in || [0J. A detailed description of this and the results for the 
corresponding numerical evaluations will be reported elsewhere. Here we just briefly outline the 
scheme with the example of the surface energy ( 3.12| ). 



First of all we note that || 

oo 

J^A = 0, D<1. (3.23) 

1=0 

As a result for these values D we can add to the integrand any /-independent term without 
changing the value of the sum. We choose this term in a way that makes the z-integral finite. 
For this we write the asymptotic behavior of the integrand in Eq.( |3.12| ) as z — > oo: 

fgW)]-^-^-^ + ..•), (3.24) 

where we introduced the notation 

o(as)_( D-3 + AA/B, B^O 

b -j D _^ B = Q . {6.2b) 

Note that the displayed terms are independent of /. Hence, without changing the value of the 
surface energy for D < 1 we can write expression ( (3.12 ) as 



27Ta to* J ™ Vz 2 -m 2 a 2 { I v {z) V \ 2z Jj 

(3.26) 

In this form the integrals are well defined for a fixed I, as for large z — * oo the subintegrand 
behaves as z~ 2 . Then we analytically continue expression ( |3.26| ) to all D. To turn the sum over / 
into a convergent series we can use procedure often used in the previous literature on the Casimir 
effect (see, for instance, [0, H, §, 0, ||). First we rescale the integration variable, z — > zv. 
Then we add and subtract the leading terms in the asymptotic expansion of the subintegrand 
for I — ► oo using the uniform asymptotic expansions for the modified Bessel functions. The 
subtracted part is finite and may be numerically evaluated, and the asymptotic part can be 
expressed via Riemann or Hurwitz zeta functions. Note that in the final result divergent terms 
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in the form of poles for the latter functions may remain. Here one should take into account that 
the local cutoff dependent parts are automatically lost in this regularization procedure. The 
corresponding terms may be important in comparing the results with experiments (for instance, 
in the electromagnetic case, see, e.g., [47, |4q]), and can be extracted by expanding the local 



expectation values (|2.34|) in powers of the distance from the boundary. 

4 Vacuum outside a sphere 

To obtain the VEV for the EMT outside a sphere we consider first the scalar vacuum in the 
layer between two concentric spheres with radii a and b, a < b. The corresponding boundary 
conditions have the form 

A r + B r -^\ (f(x) = 0, r = a,b, (4.1) 

with constant coefficients A r and B r , in general, different for the inner and outer spheres. Now 
the complete set of solutions to the field equation has the form ( |2.1C ) with the replacement 



MXr) -» g u (Xa,\r) = J u {\r)Y& (Aa) - J<f\\a)Y u {\r) , A = ^uo 2 - m 2 , (4.2) 

where Y u (z) is the Neumann function, and functions with overbars are defined in analogy to 
Eq. Q ): 



F^ a \z) = A a F(z) + B a zF'(z), A a = A a - B a n/2a, B a = Bja, a = a,b. (4.3) 



The eigenfunctions ( |4.2[) satisfy the boundary condition on the sphere r = a. From the boundary 
condition on r = b one obtains that the corresponding eigenmodes are solutions to the equation 

Cf{b/a, Aa) = 4 a) (Aa)yj b) (Aft) - jj) b) (Aft) Y"j a) (Aa) = 0. (4.4) 

The coefficients (3 a are determined from the normalization condition (|2.16| ), where now the 
integration goes over the region between the spheres, a < r < b. Using the formula for integrals 
involving the product of any two cylinder functions one obtains 

Pl = AM f\ K b (b/a,\a), (4.5) 



where N{rrik) comes from the normalization integral ( p,13|) and we use notation 

Substituting the eigenfunctions into the mode sum ( |2.8D and using addition formula ( 2.14 ) 
for the expectation value of the field product one finds 

2/ i\—n/2 00 00 

<0|^(*M*')|0> = ^^^(2Z + n)Cr /2 (cos0)^n(7,, fc )T^(Va,7^), (4.6) 

nCL D 1=0 k=l 



with 7 Vj fc = Aa being the solutions to Eq, (|4.4| ) (see also the Appendix) and 



h(z) = , = Z = a u (z, zr/a)g v {z, zr'/a)e^^/- 2 +- 2 (t'-*)_ (4 _ 7) 



V z 2 + m 2 a 2 



To sum over k we will use the summation formula ( A. 11 ). The corresponding conditions are 
satisfied if r + r' + \t — t'\ < 2b. Note that this is the case in the coincidence limit for the region 
under consideration. Applying to the sum over k in Eq.(fD| formula (|A.11| ) one obtains 



(0|^M^|0) = ^f: 7 ^Cr /2 (co S ^)( r - {a)2 md - [a)2t - (4-8) 
2naS D ^{rr') n / 2 [J jW 2 ( z ) + y, j a)2 (z) 
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rMf~ ^/„\n( a ) 

y/z 2 /a?-m?(tf - t) 



2 f°° zdz K^irjz) Gr(z,zr/a)G\?>(z,zr'/a) 

COoil 



7T J ma V^ 2 - a 2 m 2 ^)( 2 ) Ki a \z)ll b \r]z) - ( V z)ll a) (z) 
where we have introduced notations 

G^\z,y) = I u (y)Ki a \z) - li a \z)K u (y), a = a,b (4.9) 

(the function with a = b will be used below) with the modified Bessel functions. Note that we 
have assumed values A a and B a for which all zeros for Eq^d^) are real and have omitted the 
residue terms ( A.14j ). In the following we will consider this case only. 

To obtain the vacuum EMT components outside a single sphere let us consider the limit 
b — > oo. In this limit the second integral on the right of Eq. ([4.8| ) tends to zero (for large b/a the 
subintegrand is proportional to e~ 2bz l a ), whereas the first one does not depend on b. It follows 
from here that the quantity 

<0|^M*')|0) = ^^(^^W) (4.10) 

zdz g u (z, zr/a)g v {z, zr'/a) ^/ z 2/ a 2 +m 2 (1 J_ t) 
Vz 2 + m 2 a 2 jj ) a)2 (z) + Y u {a)2 (z) 

is the Wightman function for the exterior region of a single sphere with radius a. To regularize 
this expression we have to subtract the corresponding part for the unbounded space, which, as 
we saw, can be presented in the form fl2.26|). Using the relation 



g v (z,zr/a)g u (z,zr' /a) , 1^ j!?\z 



J v (zr/a)J u (zr'/a) = -I £ Jj^L H ^{zr/a)H^{zr' /a) (4.11) 



Ji a)2 {z)+Y^\z) ^ ' ' ^ ' ' 2 
with Hj?\z), a = l,2, being the Hankel functions, one obtains 



1 9/ _i_ 

(o\v(xMx')\o) - (o M \v(xMx')\o M ) = -—^r^c^i™* ) ( 4 - 12 ) 

x y dzz e J » \ Za > H^(zr)H^\zr'). 



Assuming that the function Si la \z), (Sl 2a \z)) has no zeros for < argz < -k/2 (— vr/2 < 
argz < 0) we can rotate the integration contour on the right by the angle 7r/2 for a = 1 and 
by the angle —it/2 for a = 2. The integrals over (0, ima) and (0, — ima) cancel out and after 
introducing the Bessel modified functions one obtains 

00 , li a \za) K v (zr)K„(zr') , 

dz z — i—. ; — cosh 

Ri a \za) Vz 2 - m 2 

For the VEV of the field square in the outside region this leads to 

wv " reg ™r«S D £j '7™ #<»>(*) v/^-mV 1 ' 

where D\ is defined in accord with Eq.( |2.3l| ). 



V ' z 2 — m 2 (t' — t) 
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As in the interior case the vacuum EMT is diagonal and the corresponding components can 
be presented in the form 



l(a,r) 



2vra 3 



1 °° f 
/ 

1=0 Jrr 



Ki a) {z) Vz 2 -m 



2 n 2 



q = £,p,p±, r > a, (4.15) 



where the functions F$ [f(y)] are given by relations fl2.35| ), ( p.36 ) and (2.37). As for the interior 
components, the quantities ([4.15|) diverge at the sphere surface, r = a. The leading terms of the 
asymptotic expansions are determined by same formulas ( 2.44| ) and ( |2.45| ) with the replacement 
a — r — > r — a. 

When D = 1 from Eq.( 4.15| ) we obtain the expectation values for the one dimensional semi- 
infinite region r > a. Now the consideration similar to (2.46) and (2.47) yields 



— / </: 

2vr 



\jz 2 - m 2 B a z + Ag 2z(r _ a) 
\J z 2 -m 2 B a z - A a 



p = 0, D 



(4.16) 



These results can also be obtained from Eq.( 2.47 ) and ( 2.46| ) in the limit a —* 00, r + a = const 
(in these formulas r + a is the distance from the boundary and corresponds to r — a in Eq.( f4.16| )). 
For the massless scalar Eq.( [i.l6 ) is a special case of the general formula from [41] for a Robin 
single plate in an arbitrary dimension. 

In the case of a massless scalar the asymptotic behavior for expectation values ( [4.15| ) at large 
distances from the sphere can be obtained by introducing a new integration variable y = zr/a 
and expanding the subintegrand in terms of a/r. The leading contribution for the summand 
with a given I has order (a/r) 2l+2D ~ 1 (assuming that A a ^ vB a ) and the main contribution 
comes from the / = term. As a result we can see that for r a 



D 



' P 



D - 1 



P±, 



with 



and 



D 



8r 3 lnr/a 
A a (D-2)T(D-l/2)T y — 



D±l\ a D~2 



D>2. 



(4.17) 
(4.18) 
(4.19) 



A a -B a n 2 D -^ D l 2 T 2 {D/2) r 2D ~^ 

The case D = 1 simply follows from Eq.( 4.16j ). For the conformally coupled scalar the vacuum 
densities behave as l/r 2D+l . 

In analogy to the interior case the outside subtracted surface energy density can be presented 
in the form 



p(surf) 



Or) 



SUB 



5{r-a- 0) 



4£ 



2iraSDr 



1 00 r 
1=0 Jn 



zdz 



ma Vz 2 - m 2 a 2 Kl a) (z) 



F s [K u (zr/a)}, 



(4.20) 

where F s [f(y)] is defined as in Eq.( |3.11 ). Integration of this formula gives the total surface 
energy 



E. 



(surf) _ 4£ — 1 



ext. 



2ira 



1=0 



zdz 



f(a) 



m 



2 a 2 Ri a \z) 



F s [K v {z% 



(4.21) 



localized on the outer surface of the sphere. 
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Integrating the energy density ( 4.15 ), q = e, over the volume outside a sphere we obtain the 
corresponding volume energy: 



E. 



(vol) 
ext 



1 r c 

2vra 2 Dl J 
i=o Jm 



zdz 



f(a) 



(*) 



V^ 2 -m 2 a 2 K( a )(z) 



(4.22) 



with the notation ( |3.14| ), Now the total vacuum energy for the exterior region is obtained as 
the sum of the volume and surface parts, 



E, 



ext 



E. 



(vol) 



E. 



(surf) 



2-Ka ' j mn z 

1=0 J ma 



7 {a) (z) 



0\z) 



Ft[K„(z)], 



(4.23) 



where Ft [f(y)] is defined in accordance with Eq.( 3.17j ). As we see, as in the interior case, the 
dependences on the curvature coupling £ in the surface and volume energies cancel each other. 
The radial projection of the force acting per unit area of the sphere from outside is determined 

as 



I f 

Fext = -P\r=a+ = ^D+T^ E A J 



z 3 dz ll a) (z) 



, „ oma Vz' z - m t a i K)?' \z) 



i?n? f>{a). 



F^[K u (z)]. 



Using relation (3.19), for the massless scalar field one finds 



a D S D F ( 



ext 



E 



(vol) 



+ 



AD 
4£ - 1 



(surf) 



ext 



(4.24) 



(4.25) 



As in the interior case, this relation can also be obtained directly from the continuity equa- 
tion ( 2.33j ). Recall that the second summand on the right comes from the nonzero trace for 
the nonconformally coupled massless scalar. The exterior surface, volume, total energies, and 
vacuum forces acting on the sphere written in the form ( 4.21|) -( |4,24D are divergent. Here the 
corresponding scheme to extract finite results is similar to that for the interior quantities and is 
explained in section S 

We now turn to the case of a spherical shell with zero thickness. The total vacuum energy 



including the interior and exterior contributions can be obtained by summing Eqs.( 3.15 ) and 
fl4.22p , and the resulting vacuum force by summing Eqs.( 3~T§| ) and ( 4.24p : 



£<VOl) j 



Et° l \a)+E t 



(vol). 



ext 



E{a) = E [n (a) + E ext (a), F = F in + F ext . 



(4.26) 



Assuming that for the exterior and interior boundary conditions the coefficients A and B are 
the same (this corresponds to B\ = —B a in Eqs.(|2.1]) and (|4.1| )) after some transformations the 
corresponding expressions can be presented in the forms 



£(surf)( c 



4£ 



2ira 



i=o Jn 



zdz 



ma V ' Z 2 



m 2 a 2 



+ 



n 



i v {z)k v {z) 

z~V v {z)K' v (z) 



(4.27) 



for the total surface energy, 



Ly"^ f°° 



z 2 — m 2 a 2 



2(3 + [y 1 -(3 2 + z 2 ) 



I v {z)K v (z) 

zi' v {z)k' v {z) 



(4.28) 
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for the total energy, and 

^ oo „ c 



zdz 



2/3 - ei + 



z=o 



V z 2 — m 2 a 2 



zi' v {z)k' v {z) 



(4.29) 



for the resulting force per unit area acting on the spherical shell. In these formulas we use the 
notations 

f(z) = zf 3 f(z), = A/B, (4.30) 
for a given function f(z) and £i is defined in Eq.( 2.36l ). It can be easily seen that in the case 



m 



and P = D/2 — 1 formula (4.28) coincides with the expression for the Casimir energy 



of the TM modes derived in [10]. The corresponding electromagnetic force can be obtained by 
differentiating this expression with respect to a. 
Using the relation 



2/3 + {z 2 + v 2 - (3 1 



h{z)K v {z] 



zI' u (z)K' u {z) Iu{z)K v (z) ' 

the total Casimir energy (|4.28| ) can also be presented in the form 



(4.31) 



E(a) 



— Ea 



1=0 



dz\J z 2 — rn-it 



I u (z)K u (z) 



(4.32) 



For the massless field, using the trace relation ( |3.21 ) and integrating the continuity equation 
( |2.33| ) one obtains 



E 



vo1 - a D S D F - S D I drr 
'o 



D— 2 /rpi\ 

\ J i /reg- 



(4.33) 



Note that though for the massless scalar the subintegrand is a total divergence (see, Eq.([2.5|)) 
the corresponding contribution to the right hand side of Eq.( 4.33| ) is nonzero as the function 



(d/dr)(ip 2 ) TCg is discontinous at r = a. This contribution can be expressed via the total surface 



energy, like to Eq.(|3T20|) or Eq. fl425l) . 

The formulas for the cases of Dirichlet and Neumann boundary conditions can be obtained 
taking B\ = (/3 = oo) and A\ = (/3 = 1 — D/2), respectively. For example, in the case of 
Dirichlet vacuum force per unit area one has 



F(a) 



1 



2ira D+1 S D 



OO .qo 



1=0 



zdz 



t 2 n 2 



4(n + 1)£ - n + z 



(Iu(z)K v (z))' 
I v {z)K v {z) 



(4.34) 



By using the expressions for Di and So it can be seen that in the case of a minimally coupled 
massless scalar (£ = 0, m = 0) the expression ( 4.34| ) coincides with the formula derived in 
For the Dirichlet and Neumann boundary conditions from Eq.( 4.27j ) one has 



E 



(surf) 
(D) 



-E 



(surf) 
(TV) 



2vr 



zdz 



\J z 2 — m 2 



(4.35) 



These quantities are independent of the sphere radius and, if we follow the regularization pro- 
cedure developed in 0, Gil, we conclude that for a thin spherical shell the regularized surface 
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energy vanishes for Dirichlet and Neumann scalars. More generally, for Dirichlet and Neumann 
boundary conditions the regularized interior, Eq.( |3.12| ), and exterior, Eq.( 4.2l| ), surface energies 



vanish separately. To see this note that VEV (0\f(x)d' r (p(x')\0) is zero for r = o in the Dirichlet 
case and for r' = a in the Neumann case. As a result for these values the subtracted product 
(if(x)d' r ip(x')) reg coincides with the corresponding Minkowskian part (0M\f(x)d' r (p(x')\0M) (this 
can also be seen by direct evaluation, applying to the sum over I the Gegenbauer addition the- 



orem, and taking the value of the remained integral from [45]) and is zero in the limit r' — > a 
or r — ► a if we regularize the Minkowskian part to be zero. Now we lead to our conclusion if we 
recall that the surface energy is proportional to ((p(a)d r (p(a)) TC g. This is in agreement with the 
regularization procedure from j|, 10 1. First of all, taking into account that the regularized value 



for the integro-sum on the right of Eq.( 4.35 ) is zero, we can derive from Eqs.( |3.12j ) and ( 4.21 ) 
the relations 

p (surf) _ p(surf) _ _ p (surf) _ _ p (surf) o fi \ 

in(D) — m(N) ~~ ext(D) — ext(N)- \<i.OV) 

Hence, it is sufficient to consider the Dirichlet interior energy. It can be written in the form 

i7>(surf) 

DQ 777 — T\ = lim (f(a)d r ^{r))r:eg = - lim d r (0 M \<p(a)<p(r)\0 M ), (4.37) 

where we have used that for a Dirichlet scalar {0\(p(a)d r ip(r)\0) = 0. Now taking into account 
expression ( [2.27 ) for the Minkowskian Wightmann function, expanding Km_i)/2(' m { a ~ r ))> an d 



using 

for the interior surface energy one obtains 

F (su rf ) _ g-1 [(D ^ /2] (-l) k (mar T(l + D — 2k) 

^in(D) 2 D aT(D/2) ^ k\T(D/2-k)^ Q l\ ' 1 ' 

where the square brackets denote the largest integer less than or equal to its argument. By using 
the analytic continuation over D in || it has been shown that the regularized value for the sum 
over I in Eq.flQg|) is zero, and hence this is the case for the Dirichlet surface energy as well. 



The method of extracting finite results from expressions (|427|) -(^29l) is the same as that for 



the interior and exterior quantities and is explained in section ||. In particular, subtracting the 
leading term in the asymptotic expansion of the subintegrand for large z and using Eq.(3.23), 
the total Casimir energy may be presented in the form 

E = dz^fz 2 - m 2 a 2 ^- In \z 25b °- 1 Uz)KJz)\. (4.40) 

27ra to ^ ma 

Here the individual integrals in the series are convergent. This formula can be alternatively 
derived using, for instance, the Green function method H [l(| or applying to the corresponding 
mode sum the methods already used in pll|-|p^,|]l6|- |19|1 . For the case of a D = 3 massless 
scalar with Neumann boundary condition we recover the known result from . 



5 Vacuum energy density and stresses in a spherical layer 

We have seen that in the intermediate region between two concentric spheres the Wightman 
function for the scalar field can be presented in the form (4.8), where the term with the first 



integral is the corresponding unregularized function for the region outside a single sphere with 
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radius a. The regularization of this term was carried out in the previous section. The second 
integral on the right of Eq.( [4.8D will give a finite result at the coincidence limit for a < r < b 
and needs no regularization. As a result for the regularized Wightman function in the region 
between two spheres one obtains 



(tp(x)<p{x')) T 



U 1=0 v ' 



(5.1) 



zdz 



V? 



(az, zr)G$ {az, zr') cosh \J z 2 - m 2 (t' - t) 



where (<p(x)(p(x , ))^z is the corresponding function for the vacuum outside a single sphere with 



radius a given by Eq. (pLl3D , and we have introduced the notation 

K$\bz)/Ki a \az) 



tt au (az,bz) 



Rl a \az)ir(bz) - Kh 0> (bz)I^(az) 



Mi 



(5.2) 



Substituting Eq.(5.1) into Eq^^) we obtain that the vacuum EMT has the diagonal form ( 2.32 ) 
with components 

q(a,b,r) = q(a,r) + q a (a,b,r), a < r < b, q = e,p,p ± , (5.3) 

where q(a, r) are the corresponding functions for the vacuum outside a single sphere with radius 
a. In Eq.(5.3) the additional components are in the form 



la(a,b,r) 



1 



2Tir n S D 



1=0 



z A dz 



VP 



f n au (az,bz)F^ \G^\az,zr) 



777/ 



(5-4) 

where [f(y)] is defined by relations ( [2.35D , (|2.36|) , and ( 2.37 ). The quantities ([T4|) are finite 
for a < r < b and diverge at the surface r = b. 

It can be seen that for the case of two spheres the Wightman function in the intermediate 
region can also be presented in the form 



(ip(x)ip(x')) r 



COS I 



(5.5) 



1=0 



zdz 



\J Z 2 — 777- 



- : n bl/ {az, bz)G { P (bz, zr)G® (bz, zr') cosh V z 2 — m 2 (t' — t) 



with (ip(x)cp(x')Y^ g being the regularized Wightman function for the vacuum inside a single 
sphere with radius b (see Eq.( |2.28| ) with replacement a — ► 6), and 



Qb v (az, bz) 



li a \az)/ir(bz) 



K^(az)H 0) (bz) - Rr(bz)ir(az) 



f(b) 



(b). 



f( a )/ 



(5.6) 



Note that in the coincidence limit, x' = x, the second summand on the right hand side of 
Eq.( |5.1[ ) is finite on the boundary r = a, and is divergent on the boundary r = b. Similarly, the 
second summand on the right of formula (|5.5| ) is finite on the boundary r = b and diverges on 
the boundary r = a. It follows from here that if we write the regularized Wightmann function 
in the form 

MxMz')>re g = M*M*')>£ g + ( V {x)<ptf))% + AW(x, X% (5.7) 
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then in the coincidence limit the "interference" term AW(x, x') is finite for all values a < r < b. 
Using formulas fl2.28|) and (|5.1|) it can be seen that this term may be presented as 



AW(x,x') 



where 



■1 



E21 + U „n/2, Q , 
- r — 7-6, ICOS#) 



n/2, 



KnSD ~^ (rr') n l 



zdz 



V 'z 2 — rri 



:W (a6) (r,r')cosh V z 2 -m 2 (t' - t) 



vv vi' ) -rv, ,. , -j„\ , . 



li b) (bz) Kl a) (az) 



G^\az,zr)Gl 0> (bz,zr') 
K ( u a \az)ll b) (bz) - K {b \bz)li a \az) 



(b). 



I u (zr')K u (zr) 



(5.8) 



(5.9) 



On the basis of formula ( |5.5| ) the vacuum EMT components may be written in another 
equivalent form: 



q(a,b,r)=q(b,r) + q b (a,b,r), a<r<b, q = e,p,p±, 



(5.10) 



with q(a, r) being the corresponding components for the vacuum inside a single sphere with 
radius b (expressions ( p. 34 ) with replacement a — ► 6). Here the additional components are given 
by the formula 



Qb(a,b,r) 



1 



2irr n S D 



1=0 



z 3 dz 



n, \fz r 



- 7 n bl/ (az,bz)F^ G$\bz , zr 



q = £,p, p±- 

(5.11) 



This expressions are finite for all a < r < b and diverge at the inner sphere surface r = a. 
It follows from the above that if we present the vacuum EMT components in the form 

q(a, b, r) = q(a, r) + q(b, r) + Aq(a, b, r), a < r < b, (5-12) 

then the quantities (no summation over i) 

Aq(a, b, r) = q a (a, b, r) - q(b, r) = q b (a, b, r) - q(a, r) = §uAW(x, x) (5.13) 

are finite for all a < r < b, and i = 0, 1,2 correspond to q = e,p,p±, respectively. Near the 
surface r = a it is suitable to use the first equality in Eq,( [5.13| ), as for r — > a both summands are 
finite. For the same reason the second equality is suitable for calculations near the outer surface 
r = b. In particular, the additional radial vacuum pressure on the sphere with r = a, a = a,b 
due to the existence of the second sphere ("interaction" force ) can be found from Eqs.(5.4) and 
( |5.11 ), respectively. Using the relations 







G^(rz,rz) = -B r , 

they can be presented in the form 

1 °° f 
Pa (a,b,r = a) = -——^J 



rz—G^(rz,y) \ y=rz = A r , r = a,b, 



(5.14) 



zdz 



1=0 



-Sl av (az, bz) 



(5.15) 



{a* _ 4(2? _ l)£A a B a - [zV + 1(1 + D - 2)] B 2 a } , a = a, b. 



Unlike the self-action forces this quantity is finite for a < b and needs no further regularization. 
Using the Wronskian 



Kl«\az)±lW(bz) - ll a \az)±KW(bz) 



A 2 



B 2 a (z 2 a 2 + u 2 



(5.16) 
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it can be seen that 







\A a — B a (z 2 a 2 + u 2 )] VL av (az, bz) = —n a a— In 

da 



1 lj a) \az) R (b \bz) 
ll b) (bz) Ri a \az) 



(5.17) 



where n a = 1, n& = — 1. This allows us to write the expressions ( j5.15|) for the interaction forces 
per unit surface in another equivalent form: 



p a (a,b,r = a) 



n. 



2na D - 1 S D ^ l J m V^ 2 - m 2 



zdz 



9 1 

oa 



1=0 

li a \az) R {b) (bz) 



1 - 



£,iA a B a 



I {b \bz) Ri a \az) 



Al-Bl(z 2 a 2 + v 2 ) 
a = a,b, 



(5.18) 



where £1 is defined in Eq.( |2.36| ). For Dirichlet and Neumann scalars the second term in the 
square brackets is zero. 



6 Total Casimir energy in the region between spherical surfaces 



In this section we will consider the total Casimir energy in the region between two concentric 
spherical surfaces. The expression for the unregularized VEV for the energy density can easily 



be obtained from Eq.( 4.10 ) by applying the corresponding second order operator from Eq.( |2.6| ): 

r 3 



<0|Too(x)|0) 



7T 



°° °° ^Y 3 ,T ab (^ u) 



2r,2 



(6.1) 



1=0 k=i\Xk+ ma 



where we use the notation (|3.2| ). The total volume energy in the region between the spheres is 
obtained by integrating this expression over r. Using the standard integrals involving cylinder 
functions one finds 



h a<r<b ~ ^a<r<b + 



A h B b 



Jl(lu,k) 

J 2 (blu,k/a) 



A a B a 



(6.2) 



where 



^00 00 

E a <r<b = -X E ^1 E Uv > k ' 



Uu,k = ^J^k/a 2 + m 2 , (6.3) 

Z=0 k=l 

is the total Casimir energy for the region between the plates evaluated as a mode sum of the 
zero-point energies for each normal mode. As in the case for the region inside a sphere, here 
the energy obtained by integration from the energy density differs from the energy ( |6.3| ). This 
difference can be interpreted as due to the surface energy contrbution: 



E 



(surf) _ 4£ - 1 

2a 



E^E 



Jv,kT™(jv,k 



1=0 fc=i \hlk + v 2 



A b B b 



J 2 (blv,k/a) 



A a B a 



(6.4) 



To see this, recall that there is an additional contribution to the energy density operator located 
on the boundaries r = a,b and having the form fl3.7| ). It can be shown that by evaluating 
the corresponding VEV with the help of Eq. (|4.6| ) and integrating over the region between the 
spheres we obtain formula 
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To derive expressions for the corresponding subtracted quantities we use the formula 

(T^ nvi \x))suB = -(2£ - l/2)5(x,dM)n r {d r {<p{rMr')) SUB ) r , =r (6.5) 



with the subtracted Wightman function from Eq.( |5.1|) , n r = 1 for r = a + and n r = — 1 for 
r = b — 0. Using this formula for the integrated surface energy on the surfaces r = a + (outer 
surface of the inner sphere) and r = b — (inner surface of the outer sphere) one obtains 



ETHr) = n r ^ly Di r / dZ (6.6) 

x f ll\za) Fs[Ku{rz)] + n^az,bz)Fs[G£Haz,rz)]\ , 
{Kl a \za) J 

where we use the notation ( p.ll ). With the help of the alternative representation ( |5.5| ) for the 
Wightman function in the region between the spheres we obtain another form of the surface 
energy, 

p(surf), v 4g-l ^ f°° zdz 

E * (r) = nr —^ Di L 7W=^ (6 - 7) 

x l K }.^ Zb) F s [I u (rz)] + n bl/ (az,bz)F s [G^(bz,rz)}\ , r = a + 0,b-0. 

The equivalence of formulas ([0]) and ([T7]) can be seen directly using the identities 

-^-^lF u [I u (rz)] + n bu (az,bz)Fn[G^(bz,rz)} = . 1 ° J fl, [*„(**)] 

&z)F u [G^(az,rz)] (6.8) 

with u = s,v and r = a, b. By virtue of relations ( 5.14| ) the surface energies (|6.6| ) may be 
presented in the form 

^r f) W = ^ (surf Hr)-^(^-l)I^ r ^A / -^L=n ru (az,bz), (6-9) 
r = a + 0, 6-0. 

The first term on the right corresponds to the surface energy of a single sphere and the second 
one is induced by the existence of the second sphere and is finite for a < b. For Dirichlet and 
Neumann scalars these surface energies vanish. 

Integrating the vacuum energy density ( |5.3| ) (q = e) over the region between the spheres 
with the help of integration formula ( 3.1 3j ) for the corresponding volume energy, one obtains 

= --f> T , ^ (6-10) 



L J Za) F v [K u (rz)]+n au (az,bz)F v [G^(az,rz)] 



K^iza 



Another form is obtained by using the energy density in the form ( |5.10| ) (or equivalently indentity 
fl6.8| )). Replacing the expression in braces at the upper limit r = b in ( 6. 1C|) with the help of 
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identity 
the form 



E 



(vol) 
ab 



the expression for the volume energy in the region a < r < b can be presented in 

zdz 



4 vo1) (&)+^ o1) («) + ^-£a / 

Z7T 1=0 J " 



V z 2 — m 2 



n a Q au (az, bz) \ (4£ - l)A a B a + ( 1 - ^ ) [A 2 



(6.11) 



(z 2 a 2 + u 2 ) B 2 } 



Now the total Casimir energy in the region between the spheres may be obtained by summing 
the volume and surface parts, 



E ab = E^ ol) + E^ { \a + 0) + E^\b - 0) 



1=0 J " 



10 dz 

z 



li a \za) 
Ri a \za 



(6.12) 

r=b 



■F t [K l/ (rz)]+n au (az,bz)F t [G^(az,rz)} 



where the function F t [f(y)] is defined as Eq.( |3.17j ). Using identity (|6 . 8|) , in analogy to (6.11) 
(or directly summing Eqs.(|]|) and flOU)), one obtains the formula 



E, 



ab 



oo 

E- m (b) + E ext (a)-—J2Di \ 
Zir i=o Jn 



50 dz 

z 



(6.13) 



^ n a Q au (az,bz) [A 2 a - (z 2 a 2 + v 2 ) B 2 ] , n b = -l,n a = 1. 



a=a,b 



Note that the last ("interference") term on the right is finite for a < b and needs no further 
regularization. By virtue of relation ( 5.17 ) this formula can also be presented in the form 



Eab = E in (b) + E cxt (a)-—J2 D l / dzVz 2 -m 2 — In 

71 [=Q Jm z 



li a \az) Ri b \bz) 



li b) (bz) Kl a, (az) 



(6.14) 



As we might expect the mode sum energy does not depend on the curvature coupling parameter. 

Now we turn to a system of two concentric spheres with zero thickness. In addition to the 
contribution from the intermediate region a < r < b we have to include the contributions of the 
regions r < a and r > b. In these regions the vacuum densities are the same as in the case of 
a single sphere. For the vacuum force acting per unit area of the sphere with radius r = a one 
obtains 

F a (a, b) = F(a) — n a p a (a, b,r = a), a = a,b, n a = 1, n& = —1, (6.15) 

where F(a) is the vacuum force for a single sphere with radius a (see formula ( |4.29[ )) and the 
second summand on the right corresponds to the interaction between the spheres and is given 
by Eq.( 5.15| ) or Eq.(5.18). Taking into account formula fl6.13| ) we can see that the total vacuum 
energy is equal to the sum of the vacuum energies for the separate spheres plus an additional 
"interference" term, given by the third summand on the right of ( 6,13| ). Using expression (4.32) 
for a single shell Casimir energy this leads to 



E, 



(ab) 

0<r<oo 



1=0 Jm 



d 



x In \li a) (az)K® {bz) [l^ {bz)K^ (az) - 1^ (az)K® {bz) 



(6.16) 



The divergences in this formula are the same as for single spheres with radii a and b, and can 
be extracted in the way described in section 0. 
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7 Conclusion 



We have considered the Casimir effect for a massive scalar field with general curvature cou- 
pling and satisfying the Robin boundary condition on spherically symmetric boundaries in D- 
dimensional space. Both cases of a single sphere and two concentric spherical surfaces are 
investigated. All calculations are made at zero temperature and we assume that the boundary 
conditions are frequency independent. The latter means no dispersive effects are taken into 
account. The formulation of the theory at finite temperatures can be carried out by using the 
standard analyticity properties of the finite temperature Green function (see, for instance, p0[). 
Unlike most previous studies of the scalar Casimir effect here we adopt the local approach. To 
obtain the expectation values for the energy-momentum tensor we first construct the positive 
frequency Wightmann function (note that the Wightmann function is also important in consid- 
erations of the response of a particle detector at a given state of motion) . The application of the 
generalized Abel-Plana formula to the mode sum over zeros of the corresponding combinations 
of the cylinder functions allows us to extract the Minkowskian part and to present the subtracted 
part in terms of exponentially convergent integrals. In this paper we consider a flat background 
spacetime and the subtraction of the Minkowskian part gives finite results at any strictly interior 
or exterior points in the coincidence limit. The regularized expectation values for the EMT are 
obtained by applying on the subtracted part a certain second-order differential operator and 
taking the coincidence limit. These quantities diverge as the boundary is approached. Surface 
divergences are well known in quantum field theory with boundaries and are investigated near 
an arbitrary shaped smooth boundary. They lead to divergent global quantities, such as the 
total energy or vacuum forces acting on the sphere. To regularize them and to extract numerical 
results we can apply a procedure based on analytic continuation in the dimension already used 
in H|, and briefly described in section |j| 

The expectation values of the EMT for the region inside a spherical shell are given by formulas 
( |2,34| ). These expressions are finite at interior points and diverge on the sphere surface. The 
leading term in the corresponding asymptotic expansion is the same as for the flat boundary 
case and is zero for a conformally coupled scalar. The coefficients for the subleading asymptotic 
terms will depend on the boundary curvature, Robin coefficient, and mass. The global qunatities 
for the sphere interior region are considered in section |3[ Integrating the unregularized vacuum 
energy density we show that the volume energy differs from the total vacuum energy, evaluated 
as the sum of the zero-point energies for the elementary oscillators. We argue that this difference 
is due to the additional surface contribution to the vacuum energy located on the inner surface 
of the sphere. We give the expressions for the subtracted volume and surface energies, Eqs.( 3.15 ) 
and (|3.12 ), and for the force acting per unit area of the sphere from inside, Eq.( |3.18| ). In section 
[| first we consider the scalar vacuum in the spherical layer between two surfaces. The quantities 
characterizing the vacuum outside a single sphere are obtained from this case in the limit when 
the radius of the outer sphere tends to infinity. Subtracting the parts corresponding to the 
space without boundaries, for the regularized vacuum densities we derive formulas (4.15). These 
formulas differ from the ones for the interior region by the replacements I u — > K u , K v — > I u . 
We also consider the corresponding global quantities for the outer region. As in the interior 
case there is an additional surface energy, located on the outer surface of the sphere. The 
decomposition of the total outside Casimir energy into volume and surface parts is provided. 
The regularized surface energy is zero for Dirichlet and Neumann boundary conditions. In 
this section we further consider a spherical shell with zero thickness and with the same Robin 
coefficients on the inner and outer surfaces. The expressions for the total and surface energies, 
including the contributions from the inside and outside regions, are obtained (formulas (|4.28| ) 
and ( |4.27| )). The resulting force acting per unit area of the spherical shell is also derived (formula 
( [4.29| )). In section |H| we consider the expectation values for the energy-momentum tensor in the 
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region between two spherical surfaces. The general case is investigated when the constants A^, 
Bi, i = a,b, in the Robin boundary condition are different for the outer and inner spheres. 
The expression for the total energy in this region (including the surface parts) is derived. The 
corresponding regularized Wightmann function can be presented in different equivalent forms 
(5.5) or (5.7)). Note that in the last formula the "interference" term AW(x,x') 
([5.9D is finite in the coincidence limit for all values a < r < b. The surface 



(Eqs.Q5 
defined by Eqs.(|5~8| 



divergences are contained in the two first summands on the right of Eq. (|5.7| ), corresponding 
to the interior and exterior Wightmann functions of a single sphere. The components of the 
vacuum EMT may be written in the forms ([[j^), ( |5.10| ) or ( 5.12j ). The additional (interaction) 
force acting per unit area of the sphere with radius r = a, a = a,b due to the existence of the 
second sphere is given by formula ( 5.15| ) or by (|5.18| ) and is finite for a < b. In section ^ we 
consider the total Casimir energy and its decomposition into surface and volume parts. First 
of all we show that in the general Robin case the unregularized intermediate Casimir energy 
obtained by integration of the corresponding energy density differs from the vacuum energy 
evaluated as the sum of the zero-point energies for each normal mode of frequency. As in the 
single sphere case this difference may be interpreted as due to the surface energy contribution. 
The corresponding subtracted quantities can be presented in the form ( |6.9| ), where the second 
summand on the right is due to the existence of the second sphere and is finite. The analogous 
formulas for the integrated volume energy and total Casimir energy in the intermediate region 
are in the forms ( |6. 1 1| ) and ( |6.13| ) . The latter can also be presented as Eq. ( |6. 14| ) . As in the single 
sphere case the surface contributions to the vacuum energy vanish for Dirichlet and Neumann 
scalars. Then we consider a system of two concentric spheres with zero thickness. The vacuum 
force acting on each sphere can be presented as the sum of the force acting on a single sphere 
plus an additional interaction force which is due to the existence of the second sphere. The latter 
is given by formula ( [5.15 ) or ( 5.1§| ) and is finite for a < b. Using the expression for a single 
sphere vacuum energy the total Casimir energy can also be presented in the form (6.16). For a 
D = 1 massless scalar the results given in this paper are special cases of the general formulas 
derived in Ell] for parallel plate geometry in arbitrary dimensions. 
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A Appendix: Summation formula over zeros of C" 6 (A, z) 

Here we will consider the series over zeros of the function 

Cf(X,z) = J^zj^W-Fi^)^^), (A.l) 

where the quantities with overbars are defined in accordance with Eq.( [4.3[) . This type of series 
arises in calculations of the vacuum expectation values in confined regions with boundaries 
of spherical form. To obtain a summation formula for these series we will use the generalized 
Abel-Plana formula (GAPF) pl|, 52] . In this formula as functions g(z) and f(z) let us substitute 



(rjz) 



+ 



(rjz) 



Hz) 



C2 b ( v ,z)> 



f(z) 



Hz) 



H^ a \z)0 a \z) 



(A.2) 
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where for deflniteness we shall assume that r/ > 1 and the notations F®,i = a, b, are introduced 
in accordance with Eq. (|4.3|) . The sum and difference of these functions are 



g(z) - (-l) k f(z) 



M ka \\z) h{z) 
' a \z) CZ b (ri,z) 



1,2. 



(A.3) 



The condition for the generalized Abel-Plana formula written in terms of the function h(z) is as 
follows: 

\h(z)\ < £i(x)e Cl|s/| \z\ — ► oo, z = x + iy, (A.4) 

where c\ < 2(77 — 1), x &Ba0+&B b°~ 1 e\{x) — > for x — ► +oo. Let 7^ be zeros of the function 
C® b (r], z) in the right half plane. In this section we will assume values of v , A, and B for which 
all these zeros are real and simple, and the function H$~ a \z) (H^ a \z)) has no zeros in the right 
half of the upper (lower) half-plane. For real v and A a , B a , the zeros 7^ are simple. To 
see this note that the function J u {tz)Yy a \z) — Y u (tz)jl) a \z) is a cylinder function with respect 
to t. Using the standard result for indefinite integrals containing the product of any two cylinder 
functions (sec C3. f44|) it can be seen that 



t 



J u (tz)Yj a \z)-Y u (tz)j( a \z) 



dt 



TT 2 zTf(r],z) 



Z 7f,A;i 



where we have introduced the notation 

f(<02, 



On the other hand 



JT'\z) 

4 b)2 (rjz) 



-1 



! )] - Al - Bltf - u 2 ) 



d_ C a b(ri x = ^ jV) 



(A.5) 



(A.6) 



(A.7) 



Combining the last two results we deduce that for real v, A a , B a the derivative ( A.7 ) is nonzero 
and hence the zeros z = 7^ t are simple. By using this it can be seen that 



7T , 



Res 2=7i/ k g(z) = —T% (v,7u,k)Hlu,k)- 



(Ai 



Let h{z) be an analytic function for Rez > except the poles Zk (7^ 7in), Rez^ > 0, and with a 
possible branch point at z = 0. (For the physical application in this paper (see section ||) the 
corresponding function is analytic. However to keep the formula general we will consider the 
case of a meromorphic function.) By using the GAPF, in analogy with the derivation of the 



summation formula (4.13) in |32|| , it can be seen that if it satisfies condition (A.4) and 

h(ze ni ) = -h(z) + o(z- 1 ), z^O, (A.9) 

and the integral 

p.v. 



J^ 2 ( x )+Y„ ( 



(0)2, 



(A.10) 



x 



exists, then 



fvr 2 m r b 
,i™ \ T E Hlv,k)Tf(v, lu,k) + r v [h{z)\ - p.v. / 

'^ + °° I fc=l ^° Ju 



h(x)dx 



(x) + 1? 



(o)2. 



(A.ll) 



7T 



Res 



z=o 



h(z)Hl U) (7 1 z) 



R^ fa) [h(xe wi / 2 ) + hixe-™/ 2 )] dx 



(x) Rr (x)ir (nx) - Rr (^ir ( X ) 



,(a) . 



r(&) 



(&), 
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where 7^ jm < b < 

7v,m+l- Here the functional rv[/i(2:)] is defined as 
r v [h(z)) = vr^Res ImZfc=0 



Ji a) {z)Jl b) (vz) + yj a) (z)Y u {b) (nz) h{z) 



Ji a)2 {z)+Yr\z) 
Hj lb) (rjz) h(z) 
H { J a \z) C2 b ( V ,z) 



(a)2. 



-7T 



Res 



(-l) ! lm« fc <0 



feJ=l,2 



(A.12) 



In section || formula ( A. 11 ) is used to derive the regularized vacuum energy-momentum tensor 
in the intermediate region between two spherical shells. Note that the corresponding functions 
h(z) are analytic and hence 7v[7i(z)] = for them. In the case A a = A^, B a = from Eq.( |A.ll ) 
one obtains the summation formula (4.13) in [32]. 

Note that Eq.( A.ll ) may be generalized for the functions h(z) with purely imaginary poles 
±%fc 5 Vk > 0, satisfying the condition 



h(ze m ) = -h(z) +o((zTWk) X ) , z -> ±iy k , 



(A.13) 



and in the case of the existence of purely imaginary zeros ±iy u ^, y vk > for C® b (i],z). The 
corresponding formula is obtained from Eq.( A.ll ) by adding to the right hand side residue terms 
for z = iyk, iy u ,k in the form 



7T J^Res 



z=Vk 



H?%z) h(z) 



(A.14) 



and taking the principal value of the integral on the right which exists by virtue of Eq.(A.13). 



The arguments here are similar to those for the Remark after Theorem 3 in [32| 
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